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The leptonic bound states positronium and muonium are used to constrain Galileon contributions 
to the Lamb shift of muonic hydrogen. Through the application of a variety of bounds on lepton 
compositeness, it is shown that either the assumption of equating the charge radius of a particle 
with its Galileon scale radius is incompatible with experiments or the scale of Galileons must be 
M > 1.33 GeV, too large to solve the muon problem. The possibility of stronger constraints in the 
future from true muonium is discussed. 
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I. INTRODUCTION 


lence, was shown in Ref. [3T to be 


Several measurements in muon physics ]Tj-{5] have vary¬ 
ing levels of disagreement with theoretical calculations. 
This “muon problem” could be a sign of the violation of 
lepton universality from beyond standard model (BSM) 
physics. That the muon should be more susceptible to 
new physics is intuitively reasonable from an effective 
field theory point of view, given its dramatically larger 
mass compared to the electron, which can lead to ob¬ 
servables proportional to powers of mi/ A, where A is 
the cutoff scale. This is similar to the enhancement of 
weak interactions in muonic systems discussed in Ref. [6]. 
While larger effects with the r or quarks are not impossi¬ 
ble, the relative simplicity of muonic physics (little QCD 
contamination, longer lifetime) makes signals easier to 
detect. In this paper we will consider two facets of the 
muon problem: the anomalous magnetic moment of the 
muon a and the proton charge radius from muonic hy¬ 
drogen (r P ) M . 

Over the years, a number of ideas incorporating BSM 
physics have been proposed to solve the muon prob¬ 
lem EM], A recent proposal that disformal scalars 
[25] arising from Galileon theories [20 could resolve the 
rp discrepancy [27] offers the tantalizing possibility that 
gravitational effects could be at play in atomic systems. 
We devote the work of this paper to investigating the via¬ 
bility of this model in the leptonic bound states positro¬ 
nium (e + e _ ) and muonium, (g + e~). Additionally, we 
will make speculations about the possible limits that a 
measurement of true muonium (g + g~) might set. While 
this system is yet to be detected, near-term experiments 
[28] [29] are planned to possibly detect and characterize it. 
Due to its larger reduced mass, true muonium (TM) has 
novel features compared to positronium (Ps) and muo¬ 
nium (Mu) [HIIHH]. including the capacity to produce the 
strongest constraints on Galileons. 

The most general metric, formed from only g^ v and 
a scalar field <f> that respects causality and weak equiva- 
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9pv = A{4>, X)g fl „ + X)d^(j)d v (j>, (1) 

where X = ^g^vd^ducj). The first term gives rise to 
conformal scalars, whose couplings to matter are heavily 
constrained by various fifth-force experiments. For our 
discussion, only the second term, which gives rise to the 
disformal coupling, is of importance. This Lagrangian 
interaction is given by 

c dls = B ^ X) , ( 2 ) 

where Tj 1 ' is the energy-momentum tensor of all matter 
fields as given in the Jordan frame. 

In previous work [251 GEL 33], the leading disformal 
coupling in nonrelativistic systems was shown to be a 
one-loop quantum effect that gives rise in atoms to an 
energy-level perturbation 


6E = — 


Sminij 

32t r 3 M 8 



1 

r 7 



(3) 


where m,;, m 3 are the masses of the constituent particles 
of the bound state (throughout we will use the conven¬ 
tion that mi > mj ) and M is the Galileon coupling scale. 
This potential is sufficiently singular to be divergent for 
pointlike particles like e and g. In the case of muonic hy¬ 
drogen, {gH), it is possible to handwavingly cut off the 
integral because the proton is composite, and then one 
may apply a screening mechanism [341 - 155] to prevent the 
region inside this radius from contributing. For the stan¬ 
dard, noncomposite leptons using any composite scale as 
a cutoff scale is not possible. Further we will show that in 
particular the choice of the charge radius as the Galileon 
coupling scale assumed in Ref. m is incompatible with 
what is known about leptonic bound states. 

We begin in Sec. [TT] with a discussion of how radii are 
defined in quantum field theory and the complications 
that arise with associating the Galileon cutoff radius with 
any other measured quantity. In Sec. Ill we explore ways 
to bound the possible composite scale of leptons, with 
particular emphasis on the charge radius. Section [TV] is 
devoted to deriving constraints on Galileons from Lamb 
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shift measurements in leptonic systems. Finally, we con¬ 
clude in Sec [V] with remarks on the future of bounds from 
leptonic systems. 


In analog to this, we argue that an appropriate definition 
for r, should be via a disformal form factor, and therefore 
would be 


II. ON RADII 


= -6 


dGdis 

dq 2 



( 6 ) 


In Ref. [27 , the cutoff radius, r*, was assumed to 
be approximated by the charge radius of the proton 
y/(r 2 ) = rp which is derived from hydrogen spectroscopy 
to be r P = 0.8758(77) fm. For this value of r i; M = 320 
MeV was found to explain the rp discrepancy. This 
low value for the coupling scale can be reconciled with 
other constraints by embedding the disformal scalars in 
a Galileon theory with screening mechanisms caused by 
higher-order operators[25]. With this in mind, we will 
only consider models of disformal scalars that are em¬ 
bedded in Galileon models. In this paper, we study how 
limits on lepton compositeness affect the Galileon expla¬ 
nation for the muon problem. 

There are two ways to interpret the assumption r l = 
rp. Reference m seems to use the charge radius as a 
proxy for the size of the nuclei that can be used as a 
cutoff. This interpretation is troubling because it pre¬ 
vents sensible predictions of the Galileon in other bound 
states. Since nuclei can contain neutrons with a negative 
(r 2 ) = —0.1149(35) fm 2 , the charge radius of a nuclei 
can be reduced while the size of the nuclei could increase. 
Another troubling part of using the charge radius is that 
this choice is arbitrary. A charged particle has a num¬ 
ber of radii, each one reflecting a different distribution 
(e.g. electric charge, weak interaction, neutron density, 
strange quark density, matter density), and as we will see, 
assuming any two are (nearly) equal has implications for 
other particles. A further complication of this view is 
that it provides no explanation as to how leptonic bound 
systems can regulate the divergence of Eq. ©• 

Another way to understand this assumption is that it 
expresses a relationship between the underlying distribu¬ 
tions. In this paper we investigate whether this relation¬ 
ship can be sustained quantitatively with leptonic bound 
states. To begin, formally the charge radius of a particle 
is defined via the electric form factor, 


G E (q 2 ) = J d 3 xe l ^p(x) 

= J d 3 ,( 1 + i q. a:+ ( ^ ) 2 + ...)p(x ) 

= Qtot - \\q\ 2 (r 2 ) + ■ ■■ , (4) 
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where G E is the electric form factor, p(x) is the charge 
density, and Qtot is the total charge of the particle. The 
standard definition of (r 2 ) is then 


2 > = -6 


dGt 
dq 2 


(5) 


q 2 — 0 


By this definition, we see that ri is related to a disformal 
density Pdis^) that represents the spatial distribution 
of matter coupling to the Galileons. Therefore the as¬ 
sumption Ti = rp is not an arbitrary statement, but is 
tantamount to saying 

f d 3 x\x\ 2 p(x) = J d 3 x\x\ 2 p dis (x). (7) 

Since setting the right-hand side of Eq. |7| for leptons 
to zero is unacceptable due to the divergence in Eq. (|3|, 
this implies that leptons must have a charge distribution 
different from a point particle. If instead is a prop¬ 
erty of particles unconnected to their charge radius, then 
the so-far unobserved lepton charge radius would give no 
constraint. This would be in analogy to how the Zemach 
radius is a property of charged particles arising from the 
magnetic field distribution, and therefore has no neces¬ 
sary relation to the charge radius. 


III. COMPOSITENESS OF LEPTONS 

As stated above, the results of Ref. m relied upon 
rp > 0 in order to cut off the divergences. Using the 
measured value from hydrogen or electron scattering ex¬ 
periments, a novel correction to the Lamb shift of muonic 
hydrogen can mimic a scenario with a smaller rp. Since 
Galileons couple to all matter content equally, this in¬ 
teraction should occur in leptonic systems also. But if 
leptons are truly pointlike, the potential would lack reg¬ 
ularization and would yield unacceptably large correc¬ 
tions to the Lamb shift and Is — 2s intervals. In order 
to prevent this effect, we are forced to introduce a com¬ 
posite radiu s. Here, we investigate the strict constraint 
that yj(r 2 ) = rt (i.e., the charge radius is the composite 
radius), and a more general constraint on a composite 
radius. 

Since we must demand screening mechanisms to evade 
other bounds, it is nontrivial to construct a composite 
scale from Galileons alone; therefore, we compute several 
different limits on composite lepton radii. 


A. Spectroscopy 

Derived constraints from hyperfine splitting (hfs), 
Lamb shifts, and Is — 2s intervals are either directly on 
the (r 2 ) or on the Zemach radius, (r)( 2 ). The Zemach 
radius is approximately linearly related to \J (r 2 ), with 
a model-dependent 0( 1) coefRcient(See Ref. (37 for a 
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discussion). Composite leptons have finite-size contribu¬ 
tions similar to those of the proton [38]. For the case of an 
s-state energy level, the leading finite-size contribution is 
known, 

6E=^(Za)^ 3 (r 2 ), ( 8 ) 

where n is the principal quantum number, Z is the charge 
of the particle, and /i is the reduced mass of the system. 
At this order, p states are not affected by (r 2 ), so for the 
Lamb shift, the contribution is 

SE Lamb = ^(Za)V(r 2 ). (9) 

Furthermore for the Is — 2s interval, this contribution 
yields 

<5 E ls _ 2s = ^(Za) 4 p, 3 (r 2 ). (10) 

We can derive stronger limits from the hyperfine splitting 
(hfs), where the leading-order effect is given in Ref. TTTfl : 

6E his = -2{Za)ii{r) {2) E F , (11) 

where the Fermi energy is given by 

E f = ^(Za) 4 (l + ai) — ( — 'j m e , (12) 
3 rrii \rrij J 

where at is the anomalous magnetic moment of particle 
i. In addition to the Zemach radius, there is a higher- 
order contribution directly from (r 2 ), which is found in 

Ref. 00]: 

SE his = ^ (Za) 2 In ( Za ) p 2 (r 2 )E F . (13) 

Applying these expressions to the measured energy spec¬ 
trum of positronium and muonium states, we can set lim¬ 
its on (r)( 2 ) and (r 2 ). While the Lamb shift and Is — 2s 
constraints in muonium can be applied to either the e or 
p component, the shift to the hfs is not mass symmetric, 
so these limits apply only to muons. Table |T] is devoted 
to listing the various constraints, using the experimental 
and theoretical values found in Table |TT] We point out 
that the Zemach radius constraint in muonium gives a 
numerical bound very close to that from a M , and so in 
Fig- [1] we only label the Zemach radius. 

In Table [I] we have used a; as a stand-in for the percent 
uncertainty in a measurement of energy shifts in true 
muonium. We see that, for even a 10% measurement 
of the Lamb shift or hyperfine splitting, true muonium 
would give competitive limits on the muon’s size. 


B. Anomalous magnetic moments 

In order to solve the muon problem, considering con¬ 
straints from the anomalous magnetic moments is criti¬ 
cal. From the formalism of Brodsky and Drell 02] , one 


Atom V^i s - 2 S ( m ) 

V / F>Larnb( m ) 

v / P)hf S ( m ) 

(r) 2 (m) 

Mu 

4 x KT 15 

i x icr 14 

4 x ltr 18 

1 x 1(T 18 

Ps 

7 x 1CT 15 

8 x icr 15 

2 x itr 16 

2 x 10 -15 

TM 

4^/x x 1CT 13 

8^x x 10“ 15 

yFx 10~ 18 

x x 1(T 17 


TABLE I. Constraints on the charge radius and Zemach ra¬ 
dius for leptonic systems. For true muonium (TM), x corre¬ 
sponds to the percent precision of a future measurement. 


can use the precisely measured ai to limit possible sub¬ 
structure in the leptons. Models of composite leptons 
can generically give corrections Aa; oc mi/M *, where M* 
is the scale of new physics. As pointed out in Ref. [42] 
though, these models result in a strong fine-tuning to the 
self-energy of the form Smi oc mi /M *. A more conserva¬ 
tive, and perhaps more reasonable, estimate of ai assumes 
the existence of a chiral symmetry which causes a can¬ 
cellation of the linear term in m F leaving Acq oc (jjt) • 
If we use the current best limits from experiments com¬ 
pared to theory for the values of a; ffl 021®, 

Aa e = -9.1(8.2) x 10" 13 , 

A= 287(80) x HT 11 , (14) 

and again argue that the corresponding mass scales can 
be interpreted as limits on the radius via R ~ H/M*c , we 
find flavor-dependent limits on the composite radius of 

Re <4 x 10” 19 m, 

«1 x 10" 18 m, (15) 

where, because of the discrepancy between theory and 
experiment, the results for p are a preferred scale as op¬ 
posed to a limit. The upcoming (g — 2) M experiment 05] 
anticipates a factor of 4 improvement in the measure¬ 
ment of a /i , which could improve our limit by a factor of 
2. These limits lack a clear relation to the lepton (r 2 ) but 
are strong limits on compositeness. If we seek to explain 
the complete muon problem with Galileons, a precision 
constraint from a^ would be essential to derive. 


IV. LIMITS ON GALILEONS 


It was shown in Ref. 03 tlieLt the leading correction to 
the Lamb shift of an atomic system due to a Galileon is 




248tt 3 


miV%j 

M 8 rf 


i-U z - 


6 V «o 


(16) 


where do is the Bohr radius of the atom and Ti is the 
radius of particle i at which the divergence is cut off. To 
derive constraints from this equation, we compute the 
lcr value of <5-ELamb for each leptonic system found in Ta¬ 
bic [n] Our procedure for this is to first combine the errors 
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in quadrature, and then sum this with the observed value 


of (5-ELamb- Equating these results to Eq. (16), parts of 
the parameter space of the r,; and M below the lines in 
Fig- □ are excluded at la. We find for any fixed value 
of r*, hydrogenic systems place stronger limits than lep- 
tonic ones on the value of M. While this result at first 
seems disappointing, it is important to remember that in 
hydrogenic systems one previously assumed that = r p . 
In contrast, lepton systems can take on any values of r % , 
because no composite scale has been measured. 



FIG. 1. Allowed ( ri,M ) parameter space from various atomic 
systems. The solid lines indicate current constraints from 
Eq. (161, while the dashed line for true muonium indicates a 


future measurement with 50% precision. The parameter space 
below the lines is excluded. The horizontal line is the value of 
M to explain the rp discrepancy found in Ref. m- Shaded 
regions indicate excluded values of leptonic radii for p, e and 
hfs measurements, as explained in text and summarized in 
Table |T] To aid comparisons, we have indicated rp on the 
cc-axis. 


Using the constraints on the various radii derived in 
Sec. |III[ we see that leptonic systems require the scale of 
M to be higher than that preferred by fiH. Currently, 
positronium energy shifts alone are not well enough 
known to competitively limit the charge or Zemach ra¬ 
dius, so we have not plotted them. But combining 


positronium constraints from Eq. (16) with the a e , we 


obtain the bound M a >1.2 GeV. For muonium, all lim¬ 
its on Vi are superior bounds on M to those from the 
muonic hydrogen which require r, = rp. These bounds 
vary from > 0.67 GeV to M a > 1.33 GeV. In par¬ 

ticular, we emphasize that the constraint from the charge 
radius of the muon excludes M = 320 MeV. We therefore 
conclude that leptonic systems rule out Galileons as the 
explanation of the muon problem if the radius is related 
to the charge or Zemach radius. 

Further, we notice that a fiducial measurement of the 
Lamb shift in true muonium with only 50% precision 


would improve our limit on the scale M by a factor of 2. 
In fact, due to the weak M oc E~ 1 ^ scaling of our con¬ 
straints, a mere measurement of the existence of a Lamb 
shift in true muonium would likely provide the strongest 
limit on Galileons from leptonic physics. 

Having seen that requiring (r 2 ) em = ri for leptons 
would rule out the preferred value of M from (fiH), we 
can ask whether this assumption is necessary. It is per¬ 
haps not surprising that these two scales must be differ¬ 
ent. \J (r 2 ) em is a property of particles determined by 
their charge distribution, while should in principle be 
defined in a similar way by the field distribution that 
couples to Galileons. Since the latter couples not only 
to quarks and leptons, but photons, gluons, W’s and Z’s, 
the two distributions should differ. This is analogous to 
how the magnetic distribution of the proton results in the 
Zemach radius being different from the charge radius. Al¬ 
lowing Vi ^ rp for the proton allows the muon problem 
to be solved for any value along the (fiH) line in (ti, M) 
space not excluded by limits on M from other sources, 
and since these are currently the strong compositeness 
constraints, an unobserved lepton r, is acceptable. 


A. Limits from perturbativity 

The condition that the disformal scalars avoid other 
constraints required us to embed it in a model with 
chameleon [571 [E8] and Galileon properties. In Ref. [27] . 
the full action was given by 

s = /‘ ,4 ^( 1 ^ ~ >> 2 “ 

- V(<i>) + A(^)g, u ) , 

(17) 

where we must introduce the suppression scale A to con¬ 
trol the self-interactions of the Galileons. Within this 
model, an upper limit on A can be set by requiring per¬ 
turbative unitarity to have not been violated up to cur¬ 
rently observed energies. Reference m finds that using 
LEP data, which constraints unitarity violations up to 
200 GeV, 


P /8ttM 4 \ 4/3 

2-Kmpi V y/2m e ) 


(18) 


where (3 is an 0(1) Galileon coupling. Using this bound, 
we have a relation between A and M of 


A < 6 keV 


/ M \ 16/9 
\320 MeV/ 


(19) 


Additionally, we have assumed that the Galileon is well 
approximated by a free scalar held at least down to the 
scale of the composite radius. In order to justify this 




















5 


Atom 

Obs. 

A£ el p(MHz) 

Exp. Ref. 

AEtheory (MHz) 

Theory Ref. 

<5£(MHz) 

Mu 

Lamb 

1042(23) 

ESI 

1047.490(300) 

USES] 

-5.5(230)(0.3) 


Is - 2s 2455528941.0(9.8) 

ESI 

2455528935.4(14) 

EQESI 

5.6(98)(14) 


hfs 

4463.302765(53) 

EH 

4463.30288(55) 

ESI 

-0.000115(53)(55) 

Ps 

Lamb 

13012.42(67)(154) 

ESI 

13012.41(9) 

m 

-0.01(67)(154)(9) 


Is — 2s 

1233607216.4(32) 

ESI 

1233607222.2(6) 

m 

-5.8(34)(6) 


hfs 

203389.10(74) 

ESI 

203392.411(60) 

EMI] 

-3.31(74)(48)(6) 


hfs 

203394.2(16) (13) 

ESI 

203392.411(60) 

EMI 

1.78(160)(130)(6) 

TM 

Lamb 



1.35(5) xlO 7 

ESI 

6.8(5) xlO 6 


Is — 2s 



2.55(5)xl0 11 

ESI 

1.27(5)xl0 11 


hfs 



42330577(800) (1200) 

1 

21165288(800)(1200) 


TABLE II. Experimental and theoretical values for the necessary energy shifts in leptonic systems. For the case of true 
muonium, we have used the representative value of 50% of the theoretical values for 5E 


assumption, it was found in Ref. m that 


A > — 

Vi 


/3rrii 

mpi 


1/3 


( 20 ) 


Putting Eqs. (19) and (20) together, we can obtain ad¬ 


ditional constraints on the space of (?y, M) 

Unlike the strict limits found in the previous section, 
this constraint only applies if we demand trustworthy 
perturbative solutions from Galileons. With this caveat, 
for both muons and electrons this bound is violated if 
M = 320 MeV, again ruling out the assumption that 


V. SUMMARY AND CONCLUSIONS 

In this paper, we have shown how leptonic systems of¬ 
fer competitive constraints on the scale M of Galileons. 
From the nonobservation of a lepton charge or Zemach 
radius, we can be confident that either yj (r 2 ) ^ r, or 
the Galileon scale must be M > 1.33 GeV. This result 
would be competitive with collider and astrophysical con¬ 
straints. Going beyond the assumption that the two radii 
should be related, the Galileon model is still viable for 
solving the muon problem. 

Looking forward, the improvement of these atomic con¬ 
straints is possible. Understanding the true finite-size 
effects of Galil eons re quires going beyond the crude es¬ 
timate that \J ( r 2 ) em = Ti, and accurately determining 


their relation. There are a number of ways this could 
potentially be done. The simplest theoretical, but rela¬ 
tively difficult experimentally, approach is to note that 
since s and p states have the same dependence on M but 
a different dependence on rj, measurements of the ratio 
of the absolute energy level of each state would give a 
limit on r/ alone. Another possibility is that if one could 
construct additional Galileon observable effects in a\ or 
the hfs a similar difference in dependence might arise, al¬ 
lowing us to leverage those precision experiments further. 
One could consider improving the experimental precision, 
but since M cx r i 8E ~ 1 ^, many orders of magnitude 
improvement in experiments and theory will be required 
to improve these limits. The exception is the potential 
measurement of the spectrum in true muonium. Any 
Lamb shift measurement in true muonium will improve 
these limits, and an ambitious part-per-million-level mea¬ 
surement could completely exclude Galileons. Finally, 
we note that in positronium and true muonium, the ex¬ 
istence of an annihilation channel ll —> <f> —> ll allows for 
potentially stronger limits from energy shifts or decay- 
rate limits from processes like (p, + p~) —► e + e _ . 
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